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Resume On explore le spectre d'un peigne de Dirac pondere supporte par le 
quasicristal de Thue-Morse, et on le caracterise a un ensemble de mesure nulle 
pres, au moyen de la Conjecture de Bombieri- Taylor, pour les pics de Bragg, 
et d'une autre conjecture que Ton appelle Conjecture de Aubry-Godreche- 
Luck, pour la composante singuliere continue. La decomposition de la trans- 
formee de Fourier du peigne de Dirac pondere est obtenue dans le cadre de 
la theorie des distributions temperees. Nous montrons que I'asymptotique 
de I'arithmetique des sommes p-rarefiees de Thue-Morse (Dumont; Goldstein, 
Kelly and Speer; Grabner; Drmota and Skalba,...), precisement les fonctions 
fractales des sommes de chiffres, jouent un role fondamental dans la descrip- 
tion de la composante singuliere continue du spectre, combinees a des resultats 
classiques sur les produits de Riesz de Peyriere et de M. Queffelec. Les lois 
d'echelle dominantes des suites de mesures approximantes sont controlees sur 
une partie de la composante singuliere continue par certaines inegalites dans 
lesquelles le nombre de classes de diviseurs et le regulateur de corps quadra- 
tiques reels interviennent. 

Abstract. The spectrum of a weighted Dirac comb on the Thue-Morse quasicrys- 
tal is investigated, and characterized up to a measure zero set, by means of the 
Bombieri- Taylor conjecture, for Bragg peaks, and of another conjecture that we call 
Aubry-Godreche-Luck conjecture, for the singular continuous component. The de- 
composition of the Fourier transform of the weighted Dirac comb is obtained in terms 
of tempered distributions. We show that the asymptotic arithmetics of the p-rarefied 
sums of the Thue-Morse sequence (Dumont; Goldstein, Kelly and Speer; Grabner; 
Drmota and Skalba,...), namely the fractality of sum-of-digits functions, play a fun- 
damental role in the description of the singular continous part of the spectrum, 
combined with some classical results on Riesz products of Peyriere and M. Queffelec. 
The dominant scaling of the sequences of approximant measures on a part of the 
singular component is controlled by certain inequalities in which are involved the 
class number and the regulator of real quadratic fields. 
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1 Introduction 

The ± Prouhet-Thue-Morse sequence {r]n)nef^ is defined by 

Vn = n > 0, (1.1) 

where s{n) is equal to the sum of the 2-digits no + ni + n2 + . . . in the binary 
expansion of n = ng + ni2 + n22^ + ... It can be viewed as a fixed point of the 
substitution 1^11, 1^11 on the two letter alphabet {±1}, starting with 1 
(1 stands for —1). There exists a large literature on this sequence [AMF] [Ql]. 
Let a and b be two positive real numbers such that < 6 < a. Though there 
exists an infinite number of ways of constructing a regular aperiodic point set 
of the line [Lai] from the sequence (r]„)neN, we adopt the following definition, 
which seems to be fairly canonical. Wc call Thue-Morse quasicrystal, denoted 
by Ao_5, or simply by A (without mentioning the parametrisation by a and b), 
the point set 

A := A+U(-A+) CM (1.2) 
where A^^, or simply A+, on M+, is equal to 

{0}u{/(n):= Yl {l{a + b) + ^{a-b)vm)\n = 1,2,3,...}. (1.3) 

0<m<n— 1 

The fimction / defined by (1.3) is extended to Z by symmetry: wc put 

/(O) = by convention and f{n) — —f{—n) for n G Z, n < 0. (1-4) 

For all n G Z, |/(n + 1) — /(n)| is equal either to a or 6 so that the closed 
(generic) intervals of respective lengths a and b are the two prototiles of the 
aperiodic tiling of the line M whose {f{n))nei. is the set of vertices. Though it 
is easy to check that 

(a + 6)Z c Aa,b 

and moreover that A^^f, is a Meyer set [Lai] [M] [VG], i.e. there exists a finite 
set F = {±a, ±b, ±2a, ±2b, ±a± b} such that 

Aa,b - Aa,6 C Aa,b + F, 

the Thue-Morse qiiasicrystal, and any weighted Dirac comb on it, is considered 
as a somehow mytcrious point set, intermediate between chaotic, or random, 
and periodic [AGL] [AT] [B] [GSM] [GLl] [GL2] [KIR] [Lu] [PGA] [WWVG], 
and the interest for such systems in physics is obvious from many viewpoints. 

In this note we study the spectrum of a weighted Dirac comb /i on the 
point set A^ by using arithmetic methods, more precisely by involving sum- 
of-digits fractal fonctions associated with the rarefied sums of the Thue-Morse 
sequence (Goquet [Ct], Dumont [D], Gelfond [Gd], Grabner [Grl], Newman 
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[N], Goldstein, Kelly and Specr [GKS], Drmota and Skalba [DSl] [DS2], ...). 
For this, we hold for true two conjectures which are expressed in terms of 
scaling laws of approximant measures: the Bombieri- Taylor Conjecture and 
a conjecture that we call Aubry-Godrcche-Luck Conjecture (Subsection 3.2). 
In the language of physics, the spectrum measures the extent to which the 
intensity diffracted by n is concentrated at a real number k (wave vector). It 
it can be observed by the square modulus of the Fourier transform of /x at {k} 
[Cy] [G], or eventually of its autocorrelation [H] [La2]. 

On one hand the spectrum of the symbolic dynamical system associated 
with the Prouhet-Thue-Morse sequence is known to be singular continuous: if 

2"-l 

nn{k) ■= XI "^^j exp{-2iTrjk) (1.5) 

3=0 

denotes its Fourier transform, then 

n — 1 n—1 

|^(fc)|2 = 2" (l-cos(27r2^fc)) = 2^" J[ sm'^{TT2^k) (1.6) 

is a Riesz product constructed on the sequence {2^)j>o which has the property 
that the sequence of measures 

{2-"|rM(fc)|'rffc}„>„ (1.7) 

has a unique accumulation point, its limit, for the vague topology, which is 
a singular continuous measure (Peyriere [P] §4.1, Allouche, Mendes- France 
[AMF] Appendix I, p. 337). On the other hand, Queffelec ([Q2] §6.3.2.1) has 
shown that replacing the alphabet {±1} by {0, 1} leads to a new component 
to the measure, which is discrete and exactly localizcid at the elements of 
the group ID)2 of 2-adic rational numbers in the one-dimensional torus M/27rZ, 
explicitly 

TTl 

©2:={2^— |0<m<2",n = 0,l,2,...}. (1.8) 

Therefore the spectrum of Aa^b, and more generally the Fourier transform of 
a weighted Dirac comb /i supported by Aa,6, is expected to be the sum of a 
discrete part and a singular continuous part, each of them being a function of 
a and b. In Subsection 3.1, Theorem 3.6 (the proof of which is similar to the 
proof of Theorem 4.1 in [GVG] for sets of /3- integers, with /? a quadratic 
imitary Pisot number) shows that it is the case in the context of tempered 
distributions. 

We explore the particular Dirac comb on Aa,6 for which the weights are all 
equal to 1 on the real positive line. In Subsection 3.3 we deduce the Bragg 
component of its spectrum by classical results on Riesz products of Peyriere 
[P] and Queffelec [Q2] and by the Bombieri- Taylor conjecture (Theorem 3.9). 
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In Subsection 3.4 wc show the deep relation between the p-rarefied sums of the 
Thue-Morse sequence and the singular continuous component of its spectrum. 
We use the sum-of-digits fractal functions of the rarefaction phenomenon, re- 
called in Section 4, in agreement with the Aubry-Godrcche-Luck argument 
(Subsection 3.2) to deduce, in Subsection 5.1, the singular part S of the spec- 
trum, up to a zero measure subset (Theorem 5.1). In the subsequent Sub- 
sections wc explore the extinctions properties and the questions of visibility 
of the singular component of the spectrum as a function of the sequence of 
finite approximants. This leads to characterize prime numbers p for which 
the sequences of approximant measures at the wave vectors k = ^^p; 2% have 
the property of being "size-increasing" at k] more generally, for some classes 
of prime numbers p the Cohen-Lenstra heuristics of real quadratic fields is 
required [CLl] [CL2] [CM]. In Section 6 we indicate how to deal with other 
weighted Dirac combs on the Thue-Morse quasicrystal h.a,b and prove that 
the diffraction process remains somehow invariant under the Marcinkiewicz 
equivalence relation, for the Bragg peaks, by the Bombieri- Taylor argument. 

In addition to the mathematical interest which links diffraction spectra to 
arithmetics as mentioned in the title, the present results represent the solu- 
tions of many experimental questions of physicists on Thue-Morse quasicrystals 
which remained unexplained until now. 



2 Averaging sequences of finite approximants 

Definition 2.1. An averaging sequence {Ui)i>o of finite approximants of Aa^b 
is given by a closed interval J whose interior contains the origin and a strictly 
increasing sequence {pi)i>o of positive numbers such that 

Ui = piJnAa,b, « = 0,1,2,... 

A natural averaging sequence of such approximants for the Thue-Morse 
quasicrystal is yielded by the sequence {{—Bn) U Bjv)jv>o, where 

Bn = {xeAa,b\0<x<f{N)} TV = 0,1,2,... (2.1) 

Obviously Card(Sjv) = 2N + 1 ioi N = 0,1,2, .. . 

We have the same notions for a Radon measure /i = X^nez '-^'i'^)^f(n) 
supported by Aa^;,, where lv{x) is a bounded complex-valued function, called 
weight, and where denotes the normalized Dirac measure supported by 
the singleton {/(n)}. The support of the weight fonction u) is denoted 



supp(a;) := {n e Z I u}{n) ^ 0}. 
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Denote 

AT^'^^ = U {/(")}• 

supp(aj) 

Then an averaging sequence {Ui)i>o of finite approximants of A™^''^'^-' is given 
by a closed interval J whose interior contains the origin and a strictly increasing 
sequence {pi)i>o of positive numbers such that 

Ui = p,JnAX^'^\ / = 0,1,2,.... 

Let us denote Aa,b, resp. the set of the averaging sequences of 

finite approximants of Aa,b, resp. A^"^^''"^. If the affine hull of supp(a;) is M, 
thenA,, = ^:7P('^\ 



3 Diffraction spectra 

3.1 Fourier transform of a weighted Dirac comb on Aa,b 

Let us consider the following Radon measure supported by the Thue-Morse 
quasicrystal Aa,b: 

/i = ^ a;(n)(5/(„). (3.1) 

where U!{x) is a bounded complex- valued function (weight function). Denote 
by 5 the Dirac measure, i.e. such that (5({0}) = 1. The measure (3.1) is 
translation boimded, i.e. for all compact if C M there exists ax such that 
sup^gR \ld\{K + a) < uk, and so is a tempered distribution. Its Fourier 
transform /2 is also a tempered distribution defined by 

m = nie-"'^) = ^{n)e-"'f^^^ (3.2) 

It may or may not be a measure. 

Let us now transform (3.2) using the general form of /(n), n € Z, given by 
(1.3) and (1.4). For a real number a; e M, the integer part of x is denoted by 
[x\ and its fractional part by {a;}. 

Lemma 3.1. IfnGZ is even, then f{n) = n{a + h)/2. If n>l is odd, then 
fin) = !l{a + b) + ^vn-i. (3.3) 
Proof. Obvious since, for every even integer m>0,r]m + Vm+i =0. □ 
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Wc deduce that /(n) has the general form /(n) = ain + ao + a2{x{n)} for 
all integers n > 1, where 

ai = {a + b)/2 
ao = -{a-b)/2 

a2 = 2(a — b) 

xln) = (l + ?7„_i)/4. 

Indeed, x{n) G [0, 1) is equal to its fractional part. Now, for any a; G M, 
the "fractional part" function x — > {x} is periodic of period 1 and so is the 
piecewise continuous function e~''^"2{a;} 



Jo 



g-2,r»mx^^ (3.4) 
p-ika2 _ 1 , CO (To 

= i| — = (-l)™e-''=^sinc(^fc + m7r) (3.5) 

be the coefficients of the expansion of e~*''"2{x} Fourier series. In (3.5), sine 
denotes the cardinal sine function sinc(a;) = sina;/a;. Thus 

+ CXD 

— 00 

where the convergence is punctual in the usual sense of Fourier series for 
piecewise continuous functions. Thus 

a;(n)e-*'=^(") = ^ ^(n)e-i'=[t(«+f)+i(-f)'/^-i] 

n>l n>l 

= 5^a;(n)e-^'=[t(«+^)-^]x('^c„(fc) 



^^c„(fc)a;(n)e-'=[5(«+^)- 



„■ 1+1,1-1 



>1 



Lemma 3.2. for every integer n> 1, 

s(n) = ^ — >— (mod 2). (3.7) 

Proof. Obvious. □ 

We deduce 



^ a;(n)e-^'=^(") = ^ ^ c„(yfc)a;(n)e-^'=[t(-+^)-^] 



^ gi7rm(l+s(n-l)) 



n>l mG^ri>l 
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and then, since e™ = —1, 

rneZn>l ^ ^ 



Therefore 



= ^ a;(-n)e+''=/(") + a;(0) + ^ w(n)e-^'=^(' 



n>l 



mG^ n>l 



= '^(O) + E E (^'»(^) (-^n-i)'"u;(n)e- 



+ E E {cm{-k){-rtn-irw{-n)eM^-^y-^ 



m^lj n> 1 



Let us observe that (— ryn-i)™ = 1 if m is even and that {—rjn-i)"^ = —Vn-i 
if m is odd, for all n > 1. Then 



mez 



+ E 



\n>l 



Suppose 



and denote 



{a-b)/{a + b)eQ 



9a,b ■■= gcd(a - 6, a + 6) > 1. 



(3.8) 



(3.9) 
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Lemma 3.3. If a,b G N\ {0} are such that a/b Q is an irreducible fraction, 
then ga,b = 2 if a and b are both odd, and ga,b = I if a is odd and b is even, or 
if a is even and b is odd. 

Proof. Obvious. □ 
Remark 3.4. The hypothesis (3.9) is fundamental since 



E 

n>l 



w(±n)e ^ ^ 1 



and 



( u;(±n)r/„_ie^ * 



are then periodic of period 1 and that they can be considered, under some 
assumptions (see below), as the derivative in the distributional sense of a 
function of period 1. The existence of these four new functions of period 1 
would be impossible if (3.9) would not hold. Then if (a — b)/{a + b) ^ Q 
the decomposition into a discrete part and a continuous part of the Fourier 
transform, which is deduced from them (see below), would be impossible except 
if /2(fc) is equal to 0. In other terms we can say that the condition {a — b)/{a + 
b) ^ Q makes inexistant the diffraction phenomenon by the weighted Dirac 
comb (3.1), supported by the Thue-Morse quasicrystal Aa^, this one does not 
occur whatever the weight function a; (a;) is (no concentration of intensity at 
any wave vector A:) . 

Remark 3.5. If a = 6 > and Lu{n) — 1 for n > and uj{n) — for n < —1, 
then the tempered distribution /i(fc) is an infinite sum and is reminiscent of the 
Fourier transform (1.5) of the ±1 Prouet-Thuc-Morse sequence in the usual 
way. Let us observe that it is the sections of this tempered distribution (finite 
sums obtained by truncating its "tail") which provide the Riesz products (1.6) 
when particular averaging sequences of length a power of 2 are used. We 
will approach this dependance to the choice of the averaging sequence by the 
rarefaction phenomenon (Section 4). 

Let us assume (3.9), i.e. a/b G Q. Then we make the following hypothesis: 
Hypothesis (H) we suppose that the Fourier series 



E/a;(?7.)(.: " '' ■ ''",6 \ ^/■uj{-n)c "'" \ 

\ -i7r(a+b) / a-b\ j' '^^^P' +in(a+b) ( a-b\ j' 

n>l So,!, V a+b) n>l g^^t {"' a+b) 
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resp. 

E/ Lu{n}ii„^i ( • " '' \ M-^/ Lj{-li)ll„-i( ' " '' 

I, -m(a+b) / a-b\ j' '^^Sp. 2^(^ +i7r(a+fc) / a-b\ 

n>l ga,h y a+b) n>l g^fi V* 

converges in the punctual sense to a periodic piecewise continuous function 
fu;+{x), resp. f^-{x), Un+ix), Ur,-{x), 

of period 1, which has a derivative 

fL+{x), resp. /<l,_(a;), /<I,^+(a;), fir,-{x), 

continuous bounded on the open set 

M-Up{ap,+}, resp. M - Up{ap,_}, M - Up{ap,^+}, M-Up{ap,^_}, 

where 

are the respective discontinuity points of 

/a;+(a;),resp./a;_(x), fun+{x), fuir,-{x). 

Let 

0-J9+ = /w+(ap,-i- + 0) - ^+(ap,+ - 0), 
o-p- = fu,-{ap- + 0) - fu,-{ap- - 0), 

fprj-l- = /ajr)+ (fllp,»)+ + 0) — fur]+{ap,rj+ — 0), 
^prj— ~ fu>r)—{^p,rj— ~t~ 0) f^jy—i^CLp^j^— 0), 

be the respective jumps at ap,+, resp. Op,-, ap,ri+, 0'p,ri-, of 
/u;+(a;), resp. U-{x), f^n+ix), Uri-ix). 
Then, by a classical result ([Sz], Chap. II, §2), we have the equality 

n>l ^ 



= /:+(^^)+E-.+<^^ (3-10) 

where /^_,_ has to be taken in a distributional sense. We have similar equalities 
for the three other summations. Then we are led to the following (formal) 
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formula for the Fourier transform of /x: 

mez 



i-(-ir 



pez 



mez 



l + (-l)' 



(-1)" 



/-(It' 



V 47r 



(3.11) 



within the framework of the distribution theory. 

We now consider the convergence of the part of 'jl{k) which does not contain 
the Dirac measures (with the jumps). 

Given T > and the set of Fourier coefficients {c„i(fc) | m e Z} in (3.5), 
let Jt be the space of 4-tuples {gi{k),g2{k),gz{k),gn{k)}, where each gi{k) is 
a complex valued function, such that the series 



E 

mez 



r{k) 



1 r" 1 

(5i(fc/T))-- 



{92{klT)) 



+Cm{-k) 



i+(-ir 



(93(fc/r)) 



i-(-i)' 



(54(fc/T)) 



(3.12) 



converges (in the punctual sense) to a function G(fc) which is slowly increasing 
and locally integrable in k: there exists A> Q and v > Q such that |G(A;)| < 
A|fc|'' for A: ^ +00. As a matter of fact, for any T > 0, all 4-tuples of Fourier 
exponentials e^^^ are in Jt for all A G R, by (3.6). 

If we assume that {fi+, fi_, fir,+ JLr,-} & Jt with T = Aw/ga^b then, 
assuming true the hypothesis (H), the first part of /i(fc) which does not con- 
tain the Dirac measures defines a regular tempered distribution, say lir{k). 
From (3.11), we infer that if the Fourier transform of the measure fj, can be 
interpreted as a measure too, then the first term Jir{k) may be interpreted as 
a "continuous part" of that measure. 

Summarizing we have proved the following result. 
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Theorem 3.6. Assume a/b € Q and let ga.t =gcd(a — 6, a + b). Suppose that 
the hypothesis (H) is true and that {f^+, fi-y fLri+' fLr/-} belongs to Jt with 
T = 4'K/ga,b. Then 

(i) the Fourier transfom li of /j. is the sum of a pure point tempered distri- 
bution Jlpp and a regular tempered distribution fir as follows: 



with /Ur-(fc) = w(0) 



1 — (—^ 

ifUik/T)) {fi,+ {k/T)) 



and with fippik) 



(3.13) 



1- (-1)"' 



+Cm{-k) 



2 \T J 2 

l + (-l)" fk \ l-(-l)'" 



'^pr]+°\Jf, ~ 0'P,V+ J 2 '^PV-"\'^ ~ C'P,V-) 

(3.14) 



(ii) the pure point part Jlpp is supported by the union 



3.2 Scaling behaviour of approximant measures 

It has becoming traditional for the last two decades to define the Bragg spec- 
trum, i.e. the more or less bright spots which are observed in a diffraction 
experiment on a long-range order material (say defined by Xi ^ ^ weighted 
sum of Dirac measures localized at the atomic sites) , as the pure-point com- 
ponent of the diffraction measure, defined by the Fourier transform of the 
so-called autocorrelation 7 of the measure x (for the definition of 7 in the 
context of aperiodic point sets, sec [H]). 

This calls for some assumptions: on the point set of diffractive sites (whether 
it is a model set, a Meyer set or a Dclonc set with additional properties, etc), on 
the uniqueness of the autocorrelation 7 (Theorem 3.4 in [H]), on the existence 
of a limit in the sense of Bohr-Besicovich for the averaged Fourier transform 
of finite approximants of the measure x (Theorems 5.1 and 5.3 in [H]). Under 
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these assumptions, the Bragg component of the spectrum can be computed 
using the so-called Bombieri- Taylor argument. 

Let us make precise this context in the case of a weighted Dirac comb 
H = X^nez ^a.b 9-8 (3.1). In thc general context of diffraction 

physics, the diffracted intensity at k, per diffracting site, (with a/b G Q; see 
Remark 3.4) is [Cy] [G] [H] [La2]: 



for any averaging sequence {Ui)i>o e Aa,b- 

The Bombieri-Taylor argument asserts that when (3.15) takes a nonzero 
value, for an averaging sequence (J7;);>o, then we have the equality 



where ^pp{{k}) is the value at {k} of the pure-point part of the Fourier trans- 
form 7 of the autocorrelation 7. 

The Bombieri-Taylor argument has appeared in [BTl] and [BT2] without 
proof. Since then it has become a conjecture (sec Hof [H]) since its proof 
seems fairly difficult to establish in general for Delone sets or Meyer sets [H]. 
It is only proved for model sets (Hof [H]), for some Meyer sets (Strungaru 
[Su]), and in particular contexts (de Oliveira [Oa], Lenz [Lz], . . .). We do not 
intend to prove it but to reformulate it in terms of the scaling properties of a 
family of approximant measures that we define below in a way which seems 
consistent with approaches followed by many authors (Cheng, Savit and Merlin 
[CSM], Aubry, Godreche and Luck [AGL] [Lu], Kolar, lochum and Raymond 
[KIR]). These measures will be "approximant" to the Fourier transform of 
the autocorrelation 7 but wc need not to have a precise definition of 7, or 
to know whether 7 is unique and well-defined, etc. This is a simplification 
of the problem. This reformulation presents the advantage to be employed 
in expressing the other argument, that we will call (AGL) Aubry-Godreche- 
Luck argument. The AGL argument is known to allow the computation of the 
singular continous component of the spectrum. 

The following definition is natural, see (1.7). 

Definition 3.7. Given ^ = J2nez '^^ ^o.,b and an averaging se- 

quence U = {Ui)i>o e _4*"P^('^) of finite approximants of A^"^^^'^\ we define 
the l-th approximant measure i/u,i{k)dk as 



2 




(3.15) 



I.{k) = 7pp({fc}) 



2 




dk 



(3.16) 
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Definition 3.8. Let k eR. Then 

(i) (BT) Bombieri- Taylor argument: k belongs to the set of Bragg 
peaks of the spectrum of Aa^b if and only if 

n^r ■ t nA,i{^)dk vu^i{k)dk 

< hm mf ' ' < hm sup < +oo. (3.17) 

;^+oo ca.Ta(ui)ak i^+oo card(c/;)a«; 

for at least one sequence of finite approximants U = {Ui)i>o G A^^^^^'^\ 

(ii) (AGL) Aubry-Godreche-Luck argument: k belongs to the sin- 
gular continuous part of the spectrum of Aa^b if and only if there exists 
a e (—1, 1), which depends only upon k (not upon U), such that 

< hmmt — — < hmsup t—- — — < +00 (3.18) 

1^+00 card([/i)"dA; ~ 1^+00 CBX&{Ui)°'dk ^ ' 

for all U = {Ui)i>o G A^Y'''^^ such that 

vud{k)dk 



^ liminf ■ 



*-i-cx3 card(J7;)"dfc 
or such that 



0<limsup <+oo (3.19) 



for all U = {Ui)i>o e Al"f such that 



a, 6 

m,i{k)dk 



= lim inf 



+00 ca.Td{Ui)"dk' 

The fact that the exponent a is < 1 in the (AGL)-argument, compared to 
the Bombieri- Taylor Conjecture (BT), means that the concentration of inten- 
sity in the diffraction process is "less infinite" along the singular continuous 
component than in the case of the Bragg peaks. 



3.3 The Bragg component 

In the rest of the paper, we assume a/6 G Q so that the diffraction phenomenon 
from any weighted Dirac comb on Aafi can occur, by Remark 3.4. Then 
without restriction of generality we can assume a.,h E Q (with still < 6 < a). 
Let us start by the Thue-Morse quasierystal defined by 

u){n) = if n < and uj{n) = 1 if n > 1, 

up till Section 5, and report on other weighted Dirac combs on in Section 
6. Hence supp(a;) = N \ {0}. Let 

Ui := [0, /(/)] n (A„,6 \ {0}) for I > 1, (3.20) 
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SO that {Ui)i>i e and cavd{Ui) = I. We call this sequence the canon- 

ical sequence, any other sequence of finite approximants in A^^^'^^"^ being a 
subsequence of it. 
Let fc e M and 

TOe2z TOe2z+i 

By (3.8), we have 

M(fc) = K{k){Y, e-i«»+*')"-(»-*'))'=) - K^(fc)(^ r?„_ie-*«"+'')"-("-''»'=). 

n>l n>l 

(3.21) 

Theorem 3.9. The Bragg component, denoted B^^^^, of the spectrum of the 
Thue-Morse quasicrystal given hy fjb = X^„>i '5/(n) on h.a,b is exactly the peri- 
odized normalized group of 2-adic rational numbers 



Z+^— D2 = \^^(r+^)\rGZ and0<m<2^,h = 0,l,2,...\ . 
a + b [a + b^ z"'' J 

(3.22) 



47r 
a + b 

Proof. First we write ^{k) in (3.21) as 

{K{k)+Kr,ik)){Y^ e-t((«+'')»-(»-''))'=)-2«^(fc)(^ l + Vn-l ^-.aa+b)n-(a-b))k) 
n>l n>l 



and we take the l-th approximant measures for / = 2'^ with /i = 0, 1, 2, . . . for 
k = ^q. We have 

2'- 



^,i{k)dk = ^ 



n— 1 



2 

477 



a + b 



dq. (3.23) 



We recognize the first summation as the Fourier transform of a lattice and the 
second summation as the Fourier transform of the Thue-Morse sequence on 
the alphabet {0, 1} studied by Qucffelec [Q2] §6.3.2.1 for which the spectrum 
is composed of a Bragg component and a singular continuous component (see 
(1.5) to (1.8) in the Introduction). We now use the (BT)-argument (3.17) to 
these approximant measures, allowing h to go to +00, to claim that the Bragg 
peaks arise from either the first summation or the second summation. The 
corresponding intensities can be computed from n{-^^q) and Kn{-^^q), using 
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(3.23), passing to the limit h — > +00, and may take a zero value (possible 
extinctions). □ 

3.4 Rarefied sums of the Thue-Morse sequence and 
singular continuous component of the spectrum 

From the expression (3.23) of the approximant measure with I = 2'' , from [Q2] 
§6.3.2.1 and the (AGL)-argument (3.18) (3.19), we infer that, apart from the 
Bragg component of the spectrum of the Thuc-Morse quasicrystal given by 
H — X)n>i'^/(») "-"^ ^a,bi the continuous component exists (i.e. is not trivial) 
and is only singular. 

Assume that k = {4:TTq) /{a + b) does not belong to the Bragg component 
(3.22) and that U is the canonical sequence of finite approximants (3.20) in 
the following. Then we see, when k runs over the singular continuous part of 
the spectrum, using the (AGL)-argument, that the scaling behaviour of (3.23) 
with a power of Z = 2^* is dictated and dominated by the scaling behaviour of 

the second summation J^t^i ^+''2"-^ e~^'^^"~^>^ with a power of 2''. 

More generally, when k lies in the singular continuous part of the spectrum, 
it is obvious to say that the scaling exponent a introduced in (3.18) (3.19), 
for the scaling behaviour of i/u,i{k)dk, takes only into account the scaling 
behaviour of the dominant term 

n=l 

with a (dominant) power of I, if it exists, I going to infinity, and if (i denotes 
this exponent, one has 

a = 2/3-1. 

It amounts to understand the scaling behaviour of the sums 

I 

n=l 

with a power of I, when / goes to infinity. We now precise these notions. 

Definition 3.10. LetpS N\{0, 1,2}. The subsequences ?7i , 'rii+2p, Vi+Spj--- 

are refered to as p-rarefied Thue-Morse sequences for z = 0, 1, . . . ,p — 1 {p is 
not necessarily a prime number). Their partial sums, called p-rarefied sums of 
the Thue-Morse sequence, are denoted 

Sp,i{n) = ^ rjj n=l,2,... (3.25) 

0<j<n 

j=i (modp) 



(3.24) 
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The following Lemma is obvious. 



Lemma 3.11. Let q gM. and denote Wq := {e 
subset of the unit circle \z\ = 1. Then, 



-2i7rnq 



G Z} the countable 



if a,nd only if v ^ Wq, 
#{n G Z I e-2""« =v} = { 1 if and only if v€Wq and q ^ Q, 

oo if and only if v G Wq and q G Q. 

If q e Q,v G Wq, and if we write q = ^, with t G Z,p G N \ {0}, 
gcd{t,p) = 1, then there exists an integer r G {0, 1, 2, . . . ,p — 1} such that 



{n G 



-2iwnq _ i _ 



= v} = r + pZ. 



First, let us define the following function ai{k) of G M and ^ > 1, with 
card(f/;) = I, as follows 

2 



card([/;)"'^''^ := 



1 



card(?7;) 



i-i 

3=0 



-2i7rjk 



(3.26) 



Proposition 3.12. For almost all k 

lim (k) 



(3.27) 



Proof. This is fairly classical. We have 

2 



J2 



-2i7Tjk 



j=0 



= 2^" Yi sin2(7r2^'/fc) = 2^" [| sm^Tnpj) 



j=0 



3=0 



where ipj = 2^<^0)i > 0, with ipo = k. We deduce 



a2.(fc) = 1 + — — log|sin(7r<^j)|. 
^^^^^ j=o 

By a theorem of Raikov [R] [K], for almost every initial value (po = k, we can 
change the summation into an integral and obtain the claim 



n-l 



lim (k) = 1 + lim — — - 

n— >+(X) n— ►+oonlOg2 



log|sin(7r(/5„)| 

3=0 



1 + 



log 2 



Jo 



log|sin(7ra;)|da; = — 1. 



□ 
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Corollary 3.13. The singular continuous component of the spectrum of the 
Thue-Morse quasicrystal given hy ji = X)„>i <5/(n) on Aa,6 is a subset of mea- 
sure zero ofR. 

Proof. Assume that k belongs to the singular component. Then, by the (AGL)- 
argument, there exists a G (—1, 1) such that 

< limmt , ^ .. < limsup jVrr < (3.28) 



(^+oo caiA{Ui)°'dk 



caxd{Ui)°'dk 



for all U = {Ui)i>o e A2 



such that 

„ / uu.i{k)dk 

7^ lim sup — — — — . 

l^+oo card([/i)"dfc 



Here, with the particular subsequence of finite approximants {U2r>), by (3.23) 
and Proposition 3.12, we have 



vu,2^{k)dk 

lim sup r— - — 7 — — 

„^+oo card(£/2'')"«fc 



' ' n^+oo card((72n)«dfc 







for almost all k, since lim„^oo(— <a! + Q!2''(fc)) = ~oi — 1 < 0. By the block 
structure of the Prouhet-Thue- Morse sequence we deduce that 

vu,i{k)dk 



limsup- „ 

i^+oo card{Ui)"dk 







holds for all U = {Ui)i>o G in a similar way. This is impossible, ex- 

cept possibly for the fcs which lie in a set of measure zero for which Proposition 
3.12 does not hold. □ 

Using the (AGL)-argiimcnt in a careful study of ^^pi;Q, and invoking Corol- 
lary 3.13, we will deduce in Section 5 that this singular continuous component 
is mostly 



perhaps up to a subset of R \ of measure zero. 

Now let g = - G Q, with t G Z,p e N\ {0}, and gcd(<,p) = 1. Assume that 



This implies that p 1 and that p is not a power of 2. For these values of p, 
by Lemma 3.11, we have, for any integer N > 1, 

2 



Np+1 



Np+l 
n=l 



le 



1 



Np+1 



p-i 

J2SpANp)e- 

3=0 



(3.29) 
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The equality (3.29) is the key relation for reducing the problem to the 
asymptotic behaviour of the p-rarefied sums Spj ( A^^) when goes to infinity. 
We will continue in Section 5 the characterization of the singular component 
of the spectrum, but before we need to recall some basic facts about the p- 
rarefied sums of the Thue-Morse sequence (in Section 4) when p is a prime 
number, and develop the non-prime case as in Grabner [Grl] [Gr2]. 



4 The rarefaction phenomenon 



We will follow Goldstein, Kelly and Speer [GKS] (which originates in Dumont 
[D] ) for the fractal structure of the p-rarefied sums of the Thue-Morse sequence. 
In their paper p is always an odd prime number whereas we need to consider 
the case in which p is an integer > 3 which is not necessarily prime (by (3.29)), 
and which is not a power of 2. 

In 1969 Newman [N] proved a remarkable conjecture of Moser which asserts 
that 



SsA^) > 



for any n > 1, 



namely that the 3-rarefied sequence of the Thue-Morse sequence contains a 
dominant proportion of ones. More precisely Newman proved, for all n > 1, 



3-^ 
20 



< 



nP 



< 5.3 



-13 



with /? : 



log 3 
log 4' 



(4.1) 



and that the limit lim„^+oo S^fi{n)ln^ does not exist. 

Then Coquet [Ct] improved (4.1) by the more precise statement 



53,o(n) 



log 3 

: n'°s* 



log 4 



(4.2) 



where e3_o('^) G {^,±1} and where ip^fi{x) is a continuous nowhere differen- 
tiable function of period 1 which assumes all its values in the closed interval 



[ inf , V'3,o(a:), sup ^"8,0(2;)] 
The latter is explicitly given by 



,. . , 'S'3,o(n) S'3,o(n) 
lim mi — , hm sup — 

JT, log 4 yi log 4 



log 3 

1 V°«* 2V3 55 



log 3 
log 4 



^3/ 3 ' 3 

Grabner [Grl] obtained in 1993 the Newman- type strict inequality for p = 5 
Shfi{n) > for any n > 1, 
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and a Coquet- type fractal description of S5fi{n) as 

55.o(n) = n^$5,o ( ^) + ^ for alln > 1 

\ log 16/ 5 

where ^^^{x) is also a continuous nowhere differentiable function of period 
1. We refer to [Grl] for the details. Similar results were found by Grabner, 
Herendi and Tichy for p = 17 [GHT]. 

In 1992 Goldstein, Kelly and Speer [GKS] proposed a general matrix ap- 
proach to deal simultaneously with all the sums 5pj(n), j e {0, 1, 2, . . . ,p— 1}. 
They proved that the Coquct-typc expressions arc very general, at least when 
p is an odd prime number. They obtained, for all odd prime numbers p and 
all n e N, the existence of 

• p continuous nowhere differentiable functions tppj{x),j = 0,1,2, ... ,p — 
1, of period 1, 

• an exponent P = which remarkably depends only upon p and not 
onj, 

• p error terms Ep,j (n) which arc bounded above uniformly, for some con- 
stant C > 0, for all n e N and all j £ {0,1, . . . ,p — 1}, in the sense 

C log A2 c \ ^ 1 

such that 

Sp,j{n) = n^^pj (-^^) +EpAn), (4.3) 



^ r s log 2 ^ 
where 

• r e {1, 2, 4} is an integer which remarkably depends upon p and not on 

• s is the order of 2 in the group (Z/pI^y of the invertible elements of the 
field Z/pZ (therefore depends upon p and not of j), 

• the real numbers 

Ai > A2 > . . . > 

are the moduli of the eigenvalues, in decreasing order, of the pxp matrix 

with integral coefficients. 

Proposition 4.1. Let S{n) := {Spfi{n), Sp^i{n), . . . , Sp^p-i{n)Y denote the 
vector of MP with integer entries (* for transposition). Then 

S{2'n) = MS{n). (4.4) 
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Proof. [GKS] p. 3. □ 

These authors obtain the fractal functions ijjp^j {x) by vectorial interpolation 

of (4.4) by constructing a fractal function F(a;) := {Fo{x), Fi{x), . .. , Fp_i(a;))* 
which obeys the self-similar property 

F(2"x) = M¥{x). (4.5) 

on [0, +00) and is such that F = S on N. 



Proposition 4.2. Let < 2 > denote the subgroup of {"L/pL)* generated by 2 
and let a < 2 > any of its cosets. Then the eigenvalues of M are 

€a = (-2i)^ n sin(27rj/p) (4.6) 

jeo<2> 

and satisfy 

a 

Proof Proposition 3.3 in [GKS]. □ 
As a corollary of (4.7) 

Ai > 1 Ai < 2* (4.8) 

always hold, the second largest A2 among the magnitudes of eigenvalues of M 
can be < 1 or > 1, and the first and second magnitudes of eigenvalues Ai and 
A2 can be explicitly computed for some classes of prime numbers p. These 
authors proved that the asymptotic growth of the summation Spj (n) is of the 
order of in the sense that 

-00 < liminf^^^^ < hmsup^H^ < + 00 (4.9) 

holds, for all odd prime number p and all j = 0, 1, . . . , p— 1. Let us observe that 
the p error terms Epj{n) become negligible at large n in (4.3) since /3i < /3, 
and therefore have no influence on the computation of the lower and upper 
bounds in (4.9). 

Several authors have investigated the bounds liminf and limsup in (4.9) and 
how the region between these bounds was attained by the continous fractal 
fimctions tppjix): Coquet [Ct], Dmnont [D], Drmota and Skalba [DSl] [DS2], 
Grabner [Grl], Grabncr, Herendi and Ticliy [GHT]. It appears that it is a 
non-trivial problem to decide whether the continuous function il>pj{x) has a 
zero or not. The only known examples where ippfi{x) has no zero are p = 3'^5' 
([Grl]) and p = 17 ([GHT]). Dumont [D] has shown that ^3,0(2;) and ■03,i(a;) 
have no zero, but that ^'3,2(2;) has a zero. By the same method he proved that 
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for all prime numbers p = 3 (mod 4) and j E {0, 1, . . . ,p — 1} such that the 
order of 2 in the multiplicative group (Z/pZ) is {p — l)/2, then 

ipp^j (x) has a zero. 

Let us observe [DS2] that the assertion that ippj {x) has no zero is more or less 
equivalent to the Newman-type inequality 

'^pj (^) > fo'^ almost all n (or Spj (n) < for almost all n) 

where "almost all" means "all but finitely many" . In [DS2] Drmota and Skalba 
prove (here p denotes an arbitrary integer > 3 and A'' > 1 an integer) 

p divisible by 3 or p = 4^ + 1 Spj (n) > for almost all n. 

They show that the only prime numbers p < 1000 which satisfy Spj (n) > 
for almost all n are p = 3, 5, 17,43,257,683 and give an asymptotic formula 
for their distribution at infinity. 

We now mention some results on the asymptotics of the sums Sp^i{n) when 
p > 3 is not a prime number. 

Proposition 4.3. Let 1 < ri,r2 be integers. The asymptotic behaviour of the 

sum Spfiin) with p = 3'"^5''^ is given by 

SpfiipN) = l-S3r,^o{pN) + ^Sr,^2,o{pN) + O(logiV), N>1, 

(4.10) 

where 

S^^^Ap^) = 5;r^(f^r^o(log4(piV/3)) + (|^y^Vi(ilog4(|^))+... 

. (I?) log,(|^ )) . (4.11) 

S,r.,o{pN) = ^^(p7V)/'Go(logi6(p7V/5)) + (^|^y^'Gi(ilogi6(|^))+... 

.(||)""'-''.....,^...,,|#,).-», (4,1.) 

where -Fq, -Fi, • • ■ , -PVi-i, Go, Gi, . . . , Gra-i are continuous nowhere differen- 
tiable functions of period 1, with ea^i (pA^/3''^) ,€5^2 {pN/b''"^) e {0, ±1}, and 

^ _ log3 ^_ log 5 



21og2' ' 41og2 

Proof [Grl] pp. 40-41. □ 
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Remark 4.4. Since a ~ 2\og2 > = A\og2 '"^^^^ scaling 
exponent relative to 5*3,0(3 A^), Proposition 4.3 shows that the scaHng be- 
haviour of Spfi{pN) with p = 3''i5''^ is dominated by that of 53,0(3 N), which 
is the rarefied sum of the Thue-Morse sequence relative to the prime num- 
ber 3; in this respect this prime number 3 can be qualified as "dominant" in 
the prime factor decomposition of p = 3^^ 5'^'^. Hence Proposition 4.3 tends 
to show in general that the computation of the dominant scaling exponent of 
Sp^i{pN), withp = p^iP^2 ■ ■ -Pm ' its prime factor decomposition (withpi ^ 2), 
< i < p — I, amounts to computing the scaling exponent relative to only one 
sum Sp^ ,i (n) where pj is the "dominant" prime number in the decomposition 

p = p^iP^2 ■ ■ - Pm ■ General results are still missing. 



5 The singular continous component 

5.1 Main Theorem 

Theorem 5.1. The singular continous component of the spectrum of the Thue- 
Morse quasicrystal given by ii = J2n>i ^fin) on h.a^b, is the set 

§ := \t&'L,p odd integer > 3, /i > 0, gcd{t,2^p) = 1, 

up to a subset of measure zero o/M \ ^pjQ- 

Proof. We will only consider that k G and leave apart a possible subset 

of measure zero in ]R\ ^pjQ by Corollary 3.13. There are two cases: either (i) 
^ = wit^ t e Z, J > 3 an odd integer, gcd(i,p) = 1, or (ii) = 

with t e Z, p > 3 an odd integer, h>l, gcd(t, 2^p) = 1. 

(i) Let us consider the first case. Let us assume first that p is an odd prime 
number. 

Definition 5.2. Let '■= e~^"^p. The compact subset of the complex plane 

where yj runs over the closed interval j^inf^-gjo,!] ^pji^)-: ^^Pxe[o,i] ('^^ ' 
j = 0, 1, . . . ,p — 1, is called the rarefaction domain at (the wave vector) k. 
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The rarefaction domain R(fc) is obviously facetted and may be a polygon. 
It may or may not contain the origin, or (facetted) holes. It belongs to the 
linear span of the powers of the root of unity and inherits the properties 
of the fractal fimctions ijjp,j{x), j = 0, 1, ... ,p — 1, which arc bounded. It is 
not reduced to a single point since, for all j = 0, 1, ... ,p — 1, the two bounds 
inf2;g[o,i] ippj{x) and sup3,g[o,i] i^pji^) ai'e distinct by construction. 

We now show that the (AGL)-argumcnt can be invoked, with suitable val- 
ues of a = a{k), for all such elements k. For U the canonical sequence (3.20), 
recall the expression of the l-th approximant measure, with / > 1, 



i/u,i{k)dk = y 



(,(J!L!)+,^(J!L!))(V e-^^^(--^n) 



2k*i 



n— 1 



dk. 



(5.2) 



The rarefaction domain R(A:) is related to the Z-th approximant measure i'u,i{k)dk, 
for I = Np + 1 and A'' > 1 a positive integer, as follows: by (3.29), (4.3) and 
(5.1) we have 



1 



Np+1 



Np+l 



le 



-2iTmi 



n=l 



1 



Np+1 



p-i 

Y,Sp,j{Np)e- 



Np + : 



-{Np)''\z + 0({Npf^-^) 



(5.3) 



where z G R{k), s > 1 is the smallest integer such that 2* = 1 (mod p), where 
Ai and A2 (< Ai) arc the first and second magnitudes of eigenvalues of the 
matrix M (see Section 4), where P = and where /3i = ^||| if A2 > 1, 

and /3i = if A2 < 1. These quantities are relative to the p-rarefied sums of 
the Thue-Morse sequence. Since k does not belong to the Bragg component 
^(tb second summal 
oi i'u,i{k)dk. Therefore 



B^J^l the second summation in (5.2) is dominant in the asymptotic behaviour 







i/u,Np+i{k)dk 



< liminf . 

card(L'jvp-i-i) 



'dk 



vu,Np+i{k)dk 

< hmsup — 

jv^-i-00 card^tyATp+i 



< +00 



(5.4) 



with 



(5.5) 
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and 



. r vu.Np+i{k)dk 

Imimt — — 

jv^+oo caxa{U Np+irdk 



inf 



(5.6) 



lim sup 



vu.Np+i{k)dk 
:ard(J7jVp+i)"c?fc 



sup 



>0. 



(5.7) 



It is easy to check that (5.4) holds for any other subsequence of U with the 
same exponent a = 2/3 — 1. Let us notice that 1 < Ai < 2* by (4.6), which im- 
pUes that — 1 < a < 1. Then we deduce the claim from the (AGL)-argument. 

General case: if now p > 3 is an odd integer, Cheng, Savit and Merlin 
([CSM] II. B. 4) have shown that a scaling exponent always exists. By Remark 
4.4 this scaling exponent seems to be equal, in the rarefaction phenomenon, 
to the scaling exponent of the Thue-Morse rarefied sequence relative to the 
"dominant" prime number in the prime number decomposition of p. However, 
no general result is known to the authors on this subject, except the case 
p = 3'^^5^^ with ri,r2 > 1 (i.e. Proposition 4.3 due to Grabner [Grl]). We 
now proceed as in the previous case. 

(ii) the second case can be deduced from the first one (i) as follows, now 
with / = 2", n > 1 : 



1 

2" 



2" 



E—2iTrj-r— 



2" n 



h—i , n—h—1 , 

2^ ]J sin2(7r2J-^) x 2""'' ]J sin2(7r2^-) 

j=0 ^ j=0 ^ 



h-l 



(5.8) 



We have 11^=0 i''^'^^ Wp) 7^ 0- Then, by (5.8), we can apply the asymptotic 
laws relative to the p-rarefied sums of the Thue- morse sequence as in (i). 
Therefore 

„ . . „ i'u,2"{k)dk vu,2"{k,)dk 

< limmf ' ^ ' , — — < limsup ' . — — < +00 (5.9) 

- card(?72'')"c^fc card(f/2n)"rffc ^ ' 
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with 

a = aik) = a{2^k) = = 2/3-1 = 2^^^-!, (5.10) 

^ ' ^ ' \ + bp' ^ slog2 ' ^ ' 

and 

liminf = |«,(fc)|2x2'^ n sin2(7r2^ -^)x inf \z\\ (5.11) 

limsup ''^:TJ''\'^\, =|K,(fc)|"x2^ n sin^7r2^";r^)x sup (5.12) 
„^+oS card([/2.)«dfc 1"^^' ^ 2V^ .GRC^^fe) ' ' ^ ^ 

The inequalities (5.9) hold for any other subsequence of U with the same 
exponent a = 2/3—1. The (AGL)-argument is now invoked and gives the 
claim. □ 



5.2 Extinction properties 

We say that a subsequence V = {Vi)i>o with Vi = [0,f{Ni)] n (Aa,6 \ {0}), of 
the canonical sequence U (3.20), has the extinction property at fc = ^pj^ € S 
if 

hm — — - — — = 

i^+oo card(Vi)"afc 

for the exponent a which corresponds uniquely to k in the (AGL)-argument. 

Proposition 5.3. Let k = € § with p odd. If the origin is not 

contained in the rarefaction domain R{k), then there exists no sequence of 
finite approximants V in ^^"^^^"^ which possesses the extinction property at 
k. 



Proof. Assume that such a sequence of finite approximants V = (Vi)i exists. 
Then we would have 



v\;^i{k)dk 
~+<^ caxd{Vi)°'dk 



lim inf 



= lim 



i'v,i{k)dk 
'"+00 card(yi)"dA; 



= = 



But inf^£j^(fc) \z\'^ 7^ and Kjy^^pj^^ 7^ 0. Contradiction. 



X inf \z\ 

zeR(fc) 

(5.13) 
□ 



Now let k = -^fh^ ^ ^ ^"^^ assume € R(fc). The problem of finding 
a sequence of finite approximants V = (Vi), with Vi = [0, /(iVjp)] fl {h-a,b \ 
{0}), which has the extinction property at k amounts to exhibit explicitly an 
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increasing sequence (TVj) of positive integers which satisfies 




by (4.3), and in particular sequences {Ni) for which 




as soon as i> io for a certain io. This problem concerns at the same time all 
the p functions ^pp,j (x) (whose individual cancellation properties are recalled 
in Section 4) and seems as difficult as finding particular sequences of integers 
which realize the min and the max of the fractal functions tppj (x) (see Coquet 
[Ct] forp = 3 and j = 0). 

Let us notice that most of the sequences of finite approximants have not 
the extinction property at k. 

5.3 Growth regimes of approximant measures and 
visibility in the spectrum 

Definition 5.4. Let k = -^^2^— € We say that the sequence of approxi- 



(i) size-increasing at k : if and only if a e (0, 1), 

(ii) etale at k : if and only if a = 0, 

(iii) size- decreasing at k : if and only if a e (—1, 0). 

These three regimes correspond exactly, by (5.5) and (5.10), to the three 
cases of exponents 



relative to the asymptotic laws of the p-rarefied sums of the Thue-Morse se- 
quence. 

On the other hand, it is usual in physics [Cy] [G] [Lu] to represent the 
fimction "intensity per diffracting site" Ii{k) (sec (3.15)) as a hmction of fc, 
for various values of sizes / — card(C/;) of finite approximant point sets of 
/U. In this respect, it is of common use to increase I in order to have a more 
precise understanding of the spectriun. Here, it is illusory to do so since, on 
the set of A;s at which the sequence of approximant measures is "etale" and 
"size-decreasing" , the functions 




l/2</J<l, 13=1/2, 0</3<l/2 



I Ii{k) 
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tend to zero much more rapidly than on the set of fcs at which the sequence of 
approximant measures is "size- increasing" : it suffices to compare (3.15) and 
(3.16) to observe this. This leads to a unique visibility in the spectrum, when 
I is large enough, of the singular peaks at the fcs for which the sequence of ap- 
proximant measures is "size-increasing" , and of the Bragg peaks characterized 
by the fact that limsupi^+oo7;(fc) is a nonzero constant, by the Bombieri- 
Taylor Conjecture. 

In the numerical study [PCA] most of the peaks in the spectrum, which are 
not Bragg peaks, are guessed to be labeled using "2^" at the denominator of 
the wave vector with the prime number p = 3, not with any other odd prime. 
Therefore it is important to characterize the set of odd prime numbers p for 
which the p-rarefied sums of the Thue- Morse sequence have the property that 
the sequences of approximant measures arc "size-increasing" at fc = -^ff;2^ ^ 
S. In Subsection 5.4 we ask this general question and solve it for some classes 
of prime numbers, using [GKS] . 

5.4 Classes of prime numbers 

Denote byP = {3, 5,...} the set of odd prime numbers. In the following, we 
only consider prime numbers p which are odd, so that 2 is invertible modulo 
p. The order of 2 in the multiplicative group (Z/pZ) , denoted by s, is > 1. 

Problem : Characterize the set of odd prime numbers p for which the 
exponent /3 = (log Ai)/ (slog 2) is > 1/2 (recall that Ai is the first magnitude 
of eigenvalues of the matrix M = {Sp^i-j{2^))o<i,j<p-i)- 

5.4.1 The class Vi Let Pi {p e P | s = _p - 1}. We have 

Vi = {3, 5, 11, 13, 19, 29, 37, 53, 59, 61, 67, 83, 101, 107, 131, 139, 

149, 163,173,179,181,197,...}. 

It is a conjecture of Artin [Lc] that Vi is infinite. The moduli of the eigenvalues 
of the matrix M = (5p^i_j(2P^-'^))o<i.j<p-i are 

Ai = p (which is (p — l)-fold degenerated), A2 = 0. 

The integer r is equal to 1, the error terms Epj{n) are bounded, and (4.3) 
reads as 

^'^'W - (rap) - -Z-^^- (^'"i 

with 

logp 
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with all functions il)pj{x) continuous, nowhere differentiable, of period 1. 
Proposition 5.5. The prime numbers p gPi for which (i > 1/2 are 3 and 5. 
Proof. Indeed, these values of p are given by the inequality: 

> I . (5.15) 



(p-l)log2 2 

Since x goes to zero when x tends to infinity, only a finite number 

of values of p satisfy (5.15). An easy computation gives 3 and 5. □ 

5.4.2 The class Pa.i Let ^2,1 := {p & F \ s = p=l (mod 4)}. We 
have 

7^2,1 = {17,41,97,137,193,...}. 

p — 1 

The moduli of the eigenvalues of the matrix M = (5'p,j_j(2T~))o<ij<p-i are 
Ai = e'^VP (with degeneracy > ^-^), A2 = e-'^VP e(0,l), 

where h is equal to the class number of the field Q{^/p) and e > 1 the funda- 
mental unit in the real quadratic field Q(i/p). The integer r is equal to 1, the 
error terms Epj (n) are bounded in modulus by 

2^ ^ 



and (4.3) reads as 

Sp,j{n) = n'^Vp,, ( (p!iflog2 )+^^'^^"^ ^^"^^^ 

with 

^ logp + 2/?.loge 
(p-l)log2 

with all functions ippj{x) continuous, nowhere differentiable, of period 1. 

Proposition 5.6. Only finitely many prime numbers p in 7^2, 1 give rise to the 
inequality f3 > 1/2. There is only one such prime number: p = 17. 

Proof. Let L{z,Xp) the Dirichlet L- function with character Xp (pS] p. 380). 
The discriminant of 'Q{^/p) is equal to p and we have the analytic class number 
formula ([BS] p. 385) 

2/iloge = VPi(l>Xp)- (5-17) 



30 



Jean Pierre Gazeau and Jean-Louis Verger-Gaugry 



By Hua's inequality ([H] Theorem 13.3, p. 328) 



we deduce 



Since the function 



P < 



logp+Vp(^ + l) 
(p-l)log2 

loga: + Vi(iH|£ + l) 



(5.18) 



(x-l)log2 

tends to zero when x goes to infinity, we deduce the claim. 

To deduce the possible values we notice that (5.18) intersects the line y = 
1/2 between x = 125 and 126. The Table below, obtained from [BS] pp. 
472-474, and PARl-GP, gives the first couples (p,/3) (with w = (1 + y/p)/2). 



p 


17 


41 


97 


137 


197 . . . 


h 


1 


1 


1 


1 


1 


e 


3 + 2a; 


27 + lOu 


5035 + 1138a; 


1595 + 298a; 


1637147+ 253970a; 


/3 


0.6332.. 


0.4339.. 


0.3490.. 


0.2398.. 


0.2672.. 



□ 



A complete investigation of the values of f3, for all p G V2,i, requires the 
knowledge of the class number h and the regulator log e in real quadratic fields, 
in the Cohen-Lenstra heuristics [CLl] [CL2] [CM]. 

5.4.3 The class Pa.a Let ^2,3 := {p G F \ s = p = 3 (mod 4)}. We 
have 

^2,3 = {7, 23, 47, 71, 79, 103, 167, 191, 199, . . .}. 

p— 1 

All the eigenvalues of the matrix M = (<S'p^i_j(2~2~))Q<^j<p_;^ have the same 
modulus which is equal to 

Ai = ^/p. 

The integer r is equal to 4, the error terms Epj(n) are bounded, and (4.3) 
reads as 

with 

^ logp 
(p-l)log2' 
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with all functions il)pj{x) continuous, nowhere differentiable, of period 1. 
Proposition 5.7. There exists no prime number p in V2,3 for which (3 > 1/2. 
Proof. Same proof as (5.5). □ 

5.4.4 An inequality For the other classes of prime numbers p than Vi, 

V2.1 or 1^2,3, the corresponding exponents [3 which control the rarefaction 
phenomenon have the foUowing property ([GKS] p. 14): 

^ slog2 (p-l)log2' 

It is expected that the number of such odd prime numbers for which (3 > 1/2 
is zero or eventually is finite. 



6 Other Dirac combs and Marcinkiewicz classes 

In the case of the equally weighted sum /i = X^„>]^ ^f{n) the fractahty of the 
sum-of-digits functions tppj {x) used for describing the singular continuous part 
of the spectrum of /x arises from the approximation, given by (4.5), of the so- 
lution of the matrix equation (4.4). In the general case of a weighted Dirac 
comb on the Thue-Morse quasicrystal Aa,b, the computation of scaling expo- 
nents arises from the existence of a matrix equation like (4.4) which correlates 
the block structure of the Prouhct-Thue-Morsc sequence by the powers of 2 
and the sequence of the weights (a;(n))„gz; then the computation of solutions 
by interpolation as in (4.3), with fractal powers of n, becomes possible. The 
Bragg and singular continuous components can be deduced in a similar way. 

Now let us turn to classes of weighted Dirac combs on Aa,b- Let £ be 
the space of complex-valued functions on Aa,b (or equivalently on Z through 
n f{n))- For w G £, we denote by ||«;|| the pseudo-norm ("norm 1") of 
Marcinkiewicz of w defined as 

where {Ui)i is an averaging sequence of finite approximants of Aa,b- The 
Marcinkiewicz space A4 is the quotient space of the subspace 

{gejC\ \\g\\ < +00} 

of C by the equivalence relation TZ defined by 



hn g <^ \\h-g\\ = 



(6.1) 
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(Bertrandias [Bsl] [B-VK], Vo Khac [VK]). The class of w is denoted by 
W in M. Though the definition of || • || depends upon the chosen averaging 
sequence {Ui)i the space A4 obviously does not. This equivalence relation is 
called Marcinkiewicz equivalence relation. The vector space A4 is nornied with 
\\g\\ = \\g\\, and is complete (Bertrandias [Bsl] [B-VK], Vo Khac [VK]). By 
jC°° we will mean the subspace of C of bounded weights endowed with the 
A^-topology in the sequel. 

The following proposition shows that the Bombieri- Taylor (BT) argument 
is compatible with the Marcinkiewicz equivalence relation TZ, i.e. that two 
weighted Dirac combs on Aa.b which are Marcinkiewicz-equivalent present the 
same Bragg component in the spectrum. 

Proposition 6.1. Let h,g & C such that h,gG A4. Then, for g e K, 

Ih{q) < \\hf, (6.2) 

hUg =^ h{q) = Igiq). (6.3) 

Proof. Immediate. □ 

The relation (6.3) means that the set of weighted Dirac combs on the Thue- 

Morse quasicrystal is classified by the Marcinkiewicz relation. The intensity 
function (per diffracting site) is a class function on Ai, when running over 
the Bragg component of the spectrum. 
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